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A recent numerical study by A. D. Dolgov, S. H. Hansen, S. Pastor, S. T. Pet- 
cov, G. G. RafFelt, and D. V. Semikoz (DHPPRS) [Nucl. Phys. B (in press), [Eepl 
' ph/0201287] found that complete or partial equilibrium between all active neutrino 



(N 



in 



\ flavours can be achieved before the big bang nucleosynthesis epoch via flavour os- 

' dilations, if the oscillation parameters are those inferred from the atmospheric and 



solar neutrino data, and, in some cases, if ^13 is also sizeable. As such, cosmological 

^ . constraints on the electron neutrino-antineutrino asymmetry are now applicable in 

all three neutrino sectors. In the present work, we provide an analytical treatment 

of the scenarios considered in DHPPRS, and demonstrate that their results are sta- 

Tlj" ' ble even for very large initial asymmetries. The equilibration mechanism can be 

' understood in terms of an MSW-like effect for a maximally mixed and effectively 

00 ■ monochromatic system. We also comment on DHPPRS 's choices of mixing param- 

, eters, and their handling of collisional effects, both of which could impinge on the 

Q ' extent of flavour equilibrium. 

(N ' 
O 

I. INTRODUCTION 

D , 

^ . One of the open questions in cosmology is the possibility of admitting a large relic 

neutrino-antineutrino asymmetry. For an ensemble of neutrinos and antineutrinos of flavour 
k>( ; a in thermal and chemical equilibrium, this asymmetry may be alternatively expressed in 
j_| ' terms of the chemical potential of the species. As of now, there are no direct observa- 
^ tions of the cosmic neutrino background. Thus the existence or otherwise of any sizeable ^ 

can only be established indirectly from requiring consistency with big bang nucleosynthesis 
(BBN), and from the study of the cosmic microwave background radiation (CMBR) spec- 
trum. Presently, BBN limits the chemical potential in the sector to be at most of order 1, 
while bounds on and ^i,^ are considerably less stringent |T| . A recent combined analysis of 
BBN and CMBR has generated the constraints, —0.01 < < 0.22, < 2.6, assuming 

no neutrino oscillations |]^. 

The situation changes if one considers also the effects of neutrino oscillations. In a new 
numerical study by A. D. Dolgov, S. H. Hansen, S. Pastor, S. T. Petcov, G. G. Raffelt, and 
D. V. Semikoz (DHPPRS) it was shown that, for the oscillation parameters inferred 
from the atmospheric neutrino data and the large mixing angle (LMA) solution of the solar 
neutrino problem, complete equilibrium between all active flavours is established prior to 
the onset of BBN at temperature T ~ 1 MeV.^ For other solar neutrino solutions, DHPPRS 
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^ Equilibration of the neutrino flavours for the atmospheric and LMA oscillation parameters was first 
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found that a partial equilibrium is possible if the mixing angle 613 is close to its present 
experimental limit of tan^^is ^ 0.065. Clearly, if flavour equilibrium holds, constraints on 
will apply to all three flavours. 

Central to the DHPPRS study is the inclusion of a highly nonlinear neutrino self in- 
teraction potential in the evolution equation for the ensemble. Previously, this term was 
found to give rise to synchronised vacuum oscillations characterised strongly by the initial 
conditions 0. This raises some questions: How sensitive are the DHPPRS results to the 
initial asymmetries? For instance, can the presence of a large asymmetry in the z/^ or I'-r 
sector delay its equilibration with its z/g counterpart? Why do the final asymmetries exhibit 
oscillatory behaviour in some but not all cases? Asymmetries that oscillate out of phase 
with each other are certainly not in equilibrium. 

To find the answers, we begin with an analysis of two-flavour oscillations involving in 
the solar neutrino parameter space. The knowledge gained therefrom will be applied to the 
three-flavour case, and, in particular, to investigating the role of 613. Whenever numerical 
quantities are called for, e.g., the mass squared differences, we shall adopt the values used 
in DHPPRS. These will be noted at the appropriate points. With the exception of the 
small mixing angle (SMA) solution, the solar mixing angle is always taken to be maximal in 
DHPPRS. This provides a motivation for us to examine also the more realistic case of large 
but not maximal mixing. 



II. TWO FLAVOURS 

A. Preliminary considerations 

Consider a two-state system consisting of and u^, where may be z/^, i^^, or a lin- 
ear combination thereof. We parameterise the transformation between the weak and mass 
eigenstates in vacuum with a mixing angle 6, 

f Ve \ _ ( COS 9 sin 9 \ f ^1 
J \— sin 9 cos 9 j 

where the states have masses mi^2 respectively. The same parameterisation applies to 
the z/g ^ system. 

For each momentum p, we write down the one-body reduced density matrices p {p for 
antineutrinos) and express them in terms of the function Pq (Pq) and a "polarisation" vector 
P(P): 

p = ( P""' P'A =-[Po + P-a], 

y Pxe Pxx J 2 

where P = P^jX + P^y + P^z, and a = o"a;X + cr^y + cr^z are the Pauli matrices. In this 
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notation, the and v^^ distribution functions at j9 are respectively 

= ^[P0 + Pj/cq(0), 
= ^[P0-P.]/eq(0), 



(3) 



for which we have chosen the reference distribution function /eq(0) to be of Fermi-Dirac 
form. 



1 



1 + eP/^-€ ' 



(4) 



with chemical potential ^ set to zero for all temperatures T . The distributions /p^ and /p^ 
may be similarly established from Eq. (Q) by replacing Pq and Pz with their antineutrino 
counterparts. The number density of a particle species '0 follows from taking the integral of 

over all momenta, = (l/27r^) f f^p'^dp. Note that the functions {Pq,P) and {Po,P) 
carry both momentum and time dependence unless otherwise stated. 

The evolution of P and P is governed by the quantum kinetic equations (QKEs) 



P=+ 



Am2 8V2Gpp^ 

— — B — ^ — E^^z 

2p 



Am2 8V2GfP^ 

Jj K -C/eeZ 

2p 3mf^ 



xP + y2G'F(J-J)xP, 
xP + V2Gf{J-J)xP, 



(5) 



where Am^ = — ml, B = sin2^^ x — cos2^^ z, Gp is the Fermi constant, is the 
mass of the W boson, E^e = (7/60)7r^T'^ is the electron-positron energy density, and we are 
assuming, at this stage, that the background medium does not distinguish between z/^ and 
I'r- Observe that refractive matter effects for this system are entirely CP symmetric — we 
have dropped the CP asymmetric term proportional to the difference between the charged 
lepton and antilepton number densities. This difference is expected to be of the order of 
the baryon asymmetry, and is negligible in comparison with the CP symmetric background. 
The last term in Eq. @ comprising the integrated polarisation vectors 



J 
J 



27r2 
1 

2^ 



P /eq(0) p'dp, 
P /eq(0) p'dp, 



(6) 



arises from neutrino self interaction 0, and is by definition also time-dependent. 

We have chosen to ignore the effects of collisions for now, since, for the bulk of the 
momentum distribution, the neutrino mean free path F"^ ~ {Cj^p T^)~^ is larger than an 
effective oscillation length ~ 27r(V!^2 + V^^)"^/^, where 



Am^ 

^sm2^, 
2p 



Am2 8V2Gpp 
v. = — cos 20 — K — E^ 



2p 



(7) 
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given the mass squared differences and mixing angles of tlie various solar neutrino solutions 
P). The evolution of this two-state system is therefore primarily oscillation-driven. We 
also do not consider repopulation from the background for consistency with the F = 
assumption, and hence Pq = Pq = 0. 

Note that the coUisionless approximation strictly does not apply to z/^ ^ Ur oscillations 
in the atmospheric parameter space. These oscillations become operational at higher tem- 
peratures, where frequent "non-forward" scattering on the background medium modifies 
significantly the system's dynamics (i.e., the system is "collision-driven"). One should also 
be cautious when dealing with oscillation-driven systems (such as Ug <-> u^) — residual ef- 
fects from collisions may still be sizeable, especially if there is a substantial period of large 
mixing. This issue will be revisited later in Sec. |IID| , but let us emphasise at this point 



that a full treatment of collisional effects on active-active neutrino oscillations is consider- 
ably more complicated than is implied in DHPPRS. The fact that both neutrino flavours 
can now participate in momentum-changing scattering processes gives rise to new terms in 
the evolution equation in addition to the simple damping of P^ and Py encountered in the 
active-sterile case. See, for example, Ref. for details. 

Ideally, we would like to compare the distribution functions of the two neutrino species 
concerned. When thermal and chemical equilibria prevail, these are completely specified 
by the chemical potentials = — ^p^, and = — ^p^, which are related to the neutrino- 
antineutrino asymmetries L^^ = {n^^ —np^)/n^ via 

Lu. = ^^^{^%.+ej, (8) 

where a = e,x, = 2C(3)T'^/7r^ is the photon number density, and ( is the Riemann zeta 
function. A non-thermal distribution (e.g., due to oscillations) generally has no well-defined 
chemical potential. However, if we demand only that equality between the number densities 
of Ue and (and separately, z/g and Ux) be established for flavour equilibration, then it is 
sufficient to track the evolution of the difference between L^^ and L,y^, 



= —{J.-l), (9) 
regardless of whether thermal and/or chemical equilibria are in place. 

B. Derivation of a collective evolution equation 

Let us now examine the time development of the "normalised" vector I = (J — J)/n^, for 
which we construct from Eqs. (|^) and (^ an evolution equation 



(10) 



5 



Equilibration between the two neutrino species will necessarily imply = 0. Note, however, 
that the converse is not generally true. 

Equation (^) is exact; our task is to find approximate solutions for P and P. To this 
end, we first rewrite the evolution equation in matrix form. 



d 
di 



(PA 










[( 









-v; 

K -K I + V2GFn^ 




\ h 0-4 



p. 

Py 
P. 



(V + 5)P. 

(11) 



where Vx and Vz are as defined in Eq. (|^. For the antineutrino system, the matrix V is 
replaced with —V. 

Before proceeding, observe that in the absence of the self interaction term S, the vector 
P would simply exhibit the usual "matter-suppressed" precession around the vector V = 
V^x + y^z at a rate 



(12) 



and gradually give way to vacuum oscillations as the electron-positron energy density drops 
off with the expansion of the universe. If, hypothetically, only the S term is present, P 
would precess around I at a rate 



provided that I changes sufficiently slowly with time (to be justified later). 
In the following, we shall work under the assumption that the inequality 



(13) 



(14) 



is maintained for all momenta at all times. This condition translates roughly into requiring 
that 



1 > 



~ 2 X 10-^ ^ 



y 


T 


|i| 


VMevy 



(15) 



for T > T-c, and 



1 > 



2^GfP n^\V\ 



0.12 



\/S.rn? 



T 



MeV 



(16) 



for T < Tr, where y = p/T, and 

T, ~ 19.8 y'-^ 



\/S.rn? 



MeV, 



(17) 



is the temperature at which the vacuum and electron-positron background terms become 
equal in magnitude. The quantity |I| may be regarded crudely as a measurement of the 
degree of alignment of the individual polarisation vectors P and P, such that if all neutrinos 
and antineutrinos are in fiavour eigenstates (i.e., P^^y — Px,y — 0), the magnitude of I equals 
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the difference in the Ve and neutrino-antineutrino asymmetries. It turns out that initial 
fulfilment of the requirement ([T^) tends to preserve the alignment until the self interaction 
term significantly weakens with the expansion of the universe through n^. Hence, assuming 
all P' and P', where the superscript "i" denotes initial, to be in the z direction, the reader 
may substitute in Eqs. (0) and (|16|) the relation |I| ~ |r| ~ |LJ^^ — L^yJ\i and see that 
these conditions are always met by the bulk of the momentum distribution for the Am^'s 



t 5 f^u 



lo- 



rn 



concerned, provided that the disparity between the initial asymmetries is, say 
magnitude.^ 

We solve the three coupled differential equations (|TT]) by transforming to an instantaneous 
diagonal basis defined as Q = WP, in which the evolution equation takes the form 



Q= U{V + S)U- 



uu- 



Q, 



(18) 



with W(V + S)U^^ = Diag(A;i, ^2, k^), where the eigenvalue ki is identically zero, and 
^2,3 are two imaginary numbers of equal magnitude but opposite signs. The associating 
eigenvectors qi,2,3 in terms of the original coordinates {x, y, z} constitute the columns of 
the transformation matrix 



U-^=U^ = (qi|q2|q3) 



(19) 



where qi is real, and q2,3 form a complex conjugate (c.c.) pair. In the limit tuy <^ cu^, the 
eigenvalues are 



ki = 0, A;2 = ~ iujs + O{uoy) 
and the real eigenvector qi is well approximated by 



(20) 



qi 




+ [C(u;vM)q2 + c.c] 



1+ [0{uv/uJsH + c-^] 



(21) 



where {q?, q^, qs} are the set of eigenvectors in the limit = 0. We shall not reproduce 
here the exact forms of the remaining two complex conjugate eigenvectors, but simply point 
out that q2 and qs together sweep out a plane perpendicular to qi. 

Equation (|18|) is not yet soluble; the term UIA~^ contains off-diagonal elements. However, 
these may be set to zero to facilitate calculations since the condition 



Xij 





1 1 


dl 




I 


dt 



< 1, 



(22) 



where i,j = 1,2,3 and i ^ j, is always satisfied in the cuy -C a;^ limit. This can be 
seen from Eq. (|10|), in which any one component of I is at most of order ci;v|I|; such that 
Xij ~ tuv/cu^ <^ 1 is consistent with earlier assumptions. 



^ The case of equal initial asymmetries has a priori no bearing for the purpose of the present work, although 
it could lead to some very interesting phenomena. 
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The formal solution to the now decoupled evolution equation 

~ Diag(l, e , e )Q , 

or equivalently in the original {x, y, z} basis, 



IS 



(23) 



W"'Diag(l, e 



( qi I q2 1 qs ) 



(24) 



( qi ■ P' \ 

where the superscript/subscript "i" denotes initial. Observe that the terms e^S'^sdt ^.^^ 
e~* ^ ^^'^^ lead to rapid oscillations which average to zero over the characteristic time-scale of 
qi,2.3- We retain only the time-averaged component and adopt the identity (^), whereupon 
P becomes 

P ^ (P' ■ r)i + O{uov/ujs)l + [O{uovloos)<^2 + C.C.] , (25) 
which, to the lowest order in uv/us, is a product of two quantities P' ■ I' and I carrying 



respectively the momentum and the time dependences of the original function. Equation (|25D 
pertains also to P, save for the replacement of P^ with P', and with —ujy. Substituting 
into Eq. (0), we obtain for the collective vector I a simplified evolution equation. 



I 



1 



III 277%. 



B 



2p 



(P' + P')-r/eq(0)p2rfp 



-Z,^i..(P' + P^) 



xl, 



(26) 

to the lowest order in uj^jus- 

The non-dissipative character of I's evolvement is immediately manifest in Eq. (pB]). Fur- 
thermore, the dynamics of the neutrino-antineutrino ensemble can be completely and simply 
determined from the initial conditions and from known external factors, independently of 
the evolution of the ensemble per se. Synchronised vacuum oscillations in multi-momentum 
systems subject to intense self interactions was discovered numerically in Ref. 0, and re- 
cently reinterpreted in Ref. JOT. Equation ( P^ ) is essentially a generalisation of the physical 
picture developed in Ref. [0] for pure vacuum oscillations. We opted to conduct a more 
systematic, first-principles derivation here for pedagogy. 

We take as the initial condition that all neutrinos and antineutrinos are in flavour eigen- 
states such that P^ and P' are parallel to I', and 



p r 



fi _ fi 
/eq(0) 



p' r 



fi - fi 

/eq(0) 



(27) 



up to a common sign. Two additional assumptions of thermal as well as chemical equilibria, 
that is, 



/eq (^Ve ) ' 



f- 

f- 
J l> . 



/eq( 
/eq( 



(28) 



then allow for a further approximation of Eq. (pB]) as 
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1 1 r /"Am^ 



-^ee [/eq(ej + fe^{-CJ " /eqlCJ " U'Cj] P'^P } >< I" (29) 



The integrals are now in a form opportune for the exploitation of these very useful identities: 

T3 



These, together with the conservation of the collective vector I, 

1 



III ~ ll^ ' ^ ■ 



[ifk-flJ-ifL-flJVdp 



i 3^ 



permit us to recast Eq. (^91) into a more illuminating and readily soluble form, 



2k^(e -ej + (e/-e 



"2^ 



B 



8V2G 



fP 



3m^ 



-EeeZ I X I. 



where 



2/ 



(30) 



(31) 



(32) 



(33) 



represents some "average" momentum. Equations (p2D and (|33D will form the basis of the 
discussions to follow.^ 



C. Discussions 

Equation (^) has a straightforward interpretation. Consider the terms inside the paren- 
theses. These are identically the vacuum and electron-positron background terms found in 



The method used to derive Eq. (|3^) from the QKEs is equaUy vahd had a CP asymmetric background 
been included. In fact, even if the background strongly differentiates between neutrinos and antineutrinos, 
the resulting evolution equation would still predict identical behaviours for both CP partners, contrary 
to the case when self interaction is absent. Unfortunately, this scenario is of little interest for the early 
universe. Perhaps it might find application in a supernova environment. 
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a typical single momentum evolution equation [Eq. (^) minus self interaction], and control 
the ensemble's "collective" matter-affected mixing angle, 



sin 29 r ^ 



(34) 



where Vx and Vz are evaluated for p = p. 

The characteristic momentum p defined in Eq. ([33| ) reflects, to some extent, the initial 
configuration of the ensemble. In addition to augmenting the neutrino-antineutrino asymme- 
try, a large positive chemical potential for a neutrino flavour tends to skew the distribution 
towards higher momenta, while leaving the spectrum of its less abundant antineutrino virtu- 
ally intact. A large negative chemical potential has the opposite effect. Hence the net result 
for p is that it follows more or less the trends of the more abundant species, and thereby 
grows with |^'|. For initial chemical potentials satisfying the constraint 

e/+e/^27r^ (35) 

the value of p ~ ttT originates from the mismatch between the averages (p) and {l/p)~^ 
taken over the function (P feq{C) / '^^'^ evaluated at ,^ = 0. 

The quantity exhibiting by far the strongest dependence on the initial conditions is the 
collective matter-affected oscillation length, or equivalently, the collective effective mass 
squared difference, 



Aml^ = K2p VK? + , (36) 

p=p 

where 

^ _ 3 y i.^Ue. ~ ^Vx ) ('371 

"2|7r2(e,-ej + (e/-e/)r 

For instance, Am^g vanishes for ^J,^ = — ^J,^, and oscillations are switched off completely, as 
pointed out in DHPPRS. (Note again that the present formulation is invalid for identical 
initial asymmetries.) 

An approximate solution to Eq. (^) can be established in the adiabatic limit by mapping 
I onto an instantaneous diagonal basis (in this case, isomorphic to an instantaneous mass 
basis) and setting the time derivative of the transformation matrix to zero. The resulting 
expression for the variable Iz is 

J, ~ (^2e^c29[ + s29^s2e[ cos^^ ^^^^') (^^) 
where c29 = cos 2^, s26 = sin 26', subject to the validity of the adiabatic condition 

<1. (39) 



7 = 



p=p 



at all times. For maximal mixing and the matter density profile at hand, 7 is a maximum 
at T = Tr [Eq. ( p!7D with y = y\, at which it takes on a particularly simple form: 



7l,=45° ^ 1-5 X 10-^ ( ) , (40) 
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where we have used the relation dT/dt ^ — 5.44T^/mpi, and m^i ^ 1.22 x 10^^ MeV is the 
Planck mass. 

Equation ( |38D predicts for maximal mixing an MS W- like effect, transforming from I\ 
to (plus some small amplitude oscillations) when vacuum oscillations overcome refractive 
matter effects (i.e., when c29c — c26). The temperature at which this transition takes place 
is given roughly by Tr in Eq. ( p^Tf ) with y = y. For initial chemical potentials satisfying the 
constraint (p5|), this turns out to be ~ 2.6 MeV for Am^ = 4.5 x 10"^ eV^ (LMA). Extremely 
large initial |.^|'s can lower Tr somewhat, but the connection is weak. As an illustration, the 
setting of = and C,l^ = 10 gives, for the same Am^, 1.9 MeV. In the cases of 

the LOW (Am^ = 1 x lO"'^ eV^) and the Vacuum (Am^ = 8 x 10"" eV^) solutions, the 
temperatures are ~ 0.9 MeV and ~ 0.3 MeV respectively, for reasonable |^'|'s. Evidently, 
only the LMA transition can take place well ahead of BBN. 

For the oscillation parameters of the SMA solution (Am^ = 7 x 10"^ eV^, sin 20 = 0.05), 
Iz remains close to its initial value even after the "transition" at Tj. ~ 1.9 MeV, since both 
the vacuum oscillation and matter refraction terms are predominantly in the negative z 
direction for Am^ > 0, and the usual MSW resonance condition cannot be satisfied. This 
has an important implication: in the absence of substantial coUisional damping and other 
effects not considered here, any equilibration of two neutrino species in a purely two-flavour 
scenario is, strictly speaking, an accident of maximal mixing. 

A second deciding factor on the efficacy of flavour equilibration is the adiabaticity of 
the transition from matter-suppressed to vacuum oscillations. Unlike that encountered in, 
for instance, solar neutrino analyses, the adiabaticity parameter 7 of Eq. (|39|) is strongly 
dependent on the initial conditions. Take for concreteness the case of C,l^ = 0. The reader 
may verify that for the LMA solution, the condition (|40| ) is always true if ^ 0.01, but 
can be badly violated at the transition point for the LOW Am^ unless ^ 0.1. 

In the case of maximal mixing, violation of adiabaticity at the transition point generally 
results in large amplitude "post-transition" oscillations about the equilibrium point at an 
angular frequency roughly equal to Amlg/2p. Naturally, this is quite a separate phenomenon 
from true equilibration. On the other hand, an adiabaticity parameter that evaluates to 
infinity at all times (because, for example, k — >■ 0) signifies that there is no transition at all. 
From the perspective of equilibrating two vastly different asymmetries, the requirement of 
\CuJ ~ 0-01 (assuming ^J,^ = 0) in the LMA case for a smooth transition is, by definition, 
not a major concern. We therefore dwell no further on this topic, but simply point out 
to the interested reader that there is an existing body of works devoted to non-adiabatic 



transitions in the solar interior |]lT|. The almost exponential density profile of the sun is 
of course very different from the matter density profile considered here. This will alter the 
post-transition oscillation frequency and amphtude somewhat, but not the essential physics 
of non-adiabatic effects. 



D. Large but not maximal mixing 

Previously, we saw that equilibration of z/g and in the absence of coUisional effects is 
peculiar to maximal mixing. In particular, if Am^ corresponds to that of the LMA solution 
to the solar neutrino problem, equilibrium is complete before BBN. However, the real LMA 
solution encompasses a range of mixing parameters that are merely large, but not maximal 
P). If we had used instead sin26' ~ 0.88, Eq. (|38D would predict for this scenario only a 
~ 50% reduction in i.e., a partial equilibrium. 
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The question now is: can colhsions save the scene? A full treatment of active-active 
collisions may be complicated. But for the purpose of exploring the possibilities, we may 
want to try including in the evolution equation (^) a simple active-sterile style damping 
term —D^P^ + Py) [and —D^P^ + Py) for the antineutrinos], where D = (l/2)AG|,p T^, and 
A is a positive number reflecting the relative amplitudes of the available inelastic scattering 
and elastic scattering processes that distinguish between the two flavours. 

Introduction of this "toy" damping term has a simple consequence for the collective evo- 
lution equation (P^): there is now an extra term of the form 

-5(4 + 4), (41) 

with D = (27C(3)/27r2)AG|T5 for |^^| < 1. For Amlg/2p > D, the corresponding solution 
is similar to Eq. (|SB|), but with a twist: 



Iz — e ■'^'i™*cos26'cCos26'J. + damped oscillations, (42) 
where fl ~ sin^ 29cD may be interpreted as an "instantaneous equilibration rate" [jl2 



Clearly, scales with the amount of mixing, causes the exponential in Eq. (^) to decrease 
with time, and thus contributes to equilibrating the two flavours concerned. This equili- 
brating power is especially useful if a partial equilibrium has already been reached via an 
MSW-like effect. However, the extent of the equilibrium is now dependent on the nature of 
the scattering processes. As a crude estimate, the setting of A = 1 will give exp[- ■ ■] ~ 0.16 
for sin26' ~ 0.88 at T ~ 1 MeV. Thus, together with the MSW-like transition, 4 may still 
be able to reach a value of ~ 0.08 prior to BBN. 

Before closing this section, let us stress again that it is not clear at this stage if a proper 
treatment of active-active collisions will indeed lead to an outcome similar to that afforded 
to us by the simple active-sterile picture. The latter can only serve as a rough guide. 



III. THREE FLAVOURS 

We shall not derive from first principles a three-flavour analogue of Eq. (|3^) . Such 
an exercise is perhaps not worthwhile since the coUisionless limit is strictly not applicable 
at higher temperatures. Furthermore, a complete set of QKEs for three active flavours 
incorporating all necessary collision terms has yet to be written down in a user-friendly 
form, and until we know how to handle its better established two-flavour counterpart, we 
shall not dwell on the fine details of the three-flavour case. 

However, granted that the role of the self interaction term is to cause the ensemble to 
behave in an effectively monochromatic manner, a qualitative, or even semi- quantitative, 
picture is still available if we suppose that, like the two- flavour case, the matter- affected 
mixing structure of the three-flavour system is determined by the single mode Hamiltonian 
(in flavour space). 

If/ O^ioWt-^ E,,0 , (43) 

substituted with p = p ~ ttT for small |^'|'s. Here, the quantity E^^^ = 
4m^(m^T/27r)^/^ exp(— m^/T) is the muon-antimuon energy density, with as the muon 
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mass, ml 2 3 are the squared masses of the three mass eigenstates, and the transformation 
between the weak and mass bases is parameterised with three Euler angles, 

1 \ / Ci3 Si3 \ / C12 S12 

U=\ C23 S23 ]{ 10 -su C12 I , (44) 



S23 C23 / \-Si3 Ci3 / \ 1 



where Cij = cos 6ij, and Sij = sin 6ij for ij = 12, 23, 13. Under this scheme, we identify the 
atmospheric and solar neutrino oscillation parameters as 

^"^sun = ^111-1711, 6'sun = di2, (45) 

and the third angle ^13 is subject to the constraint tan^ 613 ^ 0.065 from a combined analysis 
of the solar, atmospheric and CHOOZ data [If. 

Hamiltonian (^) is simple to decipher, thanks to the inherent mass hierarchy Am^^j^ 3> 
Arrig^jj (where Am^^^ ~ 3 x 10~^ eV^, and Amg^^^ < 10^^ eV^), and the fact that each flavour 
receives from the background medium a different contribution to their effective masses. The 
evolution of the three-flavour system follows essentially the dynamics of three separate and 
effectively two-flavour subsystems, and the present parameterisation of the transformation 
matrix U turns out to be very convenient for their description [jl3| . If all three mixing angles 
are nonzero, we have exactly three potentially equilibrating transitions: 

(i) The first occurs at T ~ 12 MeV, when Am^tm/2p — (8A/2G'FP/3m^)i?^^, and i/^ ^ I'r 
oscillations cease to be matter-suppressed. In a collisionless environment, the subsystem 
would undergo an MSW-like transformation which turns z/^ and Ur into the states Ux and 



l^y, 



Z/^ . Z/^, = -^(Z/^ - Z/^), Ur > Uy = ^{U^ + Ur), (46) 



assuming ^^atm = 45°. Collisions play the role of breaking these states into an incoherent 1:1 
mixture of and z/^. 

(ii) The second transition involves the states 

Ue > Ci3Z/e - SisUy, Vy > Si^Ve + Ci^^y, (47) 

and can only be realised for a nonzero ^13. This happens when (Am^tm + (^2^''^^VLn) / '^P — 
{8V2GFp/Sml,){Eee + at T ~ 5.2 MeV. (The factor 1/2 accompanying E^^^ comes 

from the fact that the state Uy is only "half-sensitive" to the muon-antimuon background.) 
In the absence of collisions, Eq. (^) suggests this transition to be quite impotent, even 
when ^13 is at its upper limit. However, noting all the caveats regarding active-active 
collisions, a heuristic approach to collisional damping in the style of Sec. |IID| can yield for 



tan^ 9i3 ~ 0.065 at T ~ 1 MeV an 1^ that is less than 1% of its original value by Eq. (^ 
Thus, near equilibrium between L^^ and L^^^ prior to BBN is very probable, 
(iii) The third transition, 

Cisl'e - SnUy > Ci2(Ci3Z/e - Si^Uy) - SuU^, 

^ Si2(Ci3Z/e - Si'iVy) + Ci2Z/^, (48) 
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was described in detail in Sec. || for ^13 = 0. In the case of a sizeable ^13, equilibrium between 
and L^^ should already be partially accomplished by this time. A second equilibrating 
transition between i/g and at this point will bring L^^, L^^ and L^^ even more in line. 

Note that we have been very specific with the labelling of and Vy. This actually has an 
interesting consequence. Consider the case of a maximal 6'sun and a vanishing ^13. We may 
attempt to use a simple "counting" method to establish crudely the final asymmetries given 
a set of initial conditions. For example, for D^^ = L\,^ = 0, and L]^^ = 0.1, the first transition 
distributes the asymmetry in equally amongst z/^ and Ur such that L^^ ^ Lj,^ ~ 0.05. 
The second transition is not present. At the third transition, the asymmetry carried by 

= {^fi — ^t)/V^ is shared evenly with z/g, leading to L,^^ ~ L^^ ~ 0.025. However, the 
decoupled state Uy = (z/^ + z/^)/-\/2 still has an asymmetry of 0.05. Thus the real L^,^ and L^^ 
should be ~ 0.0375. This simple counting exercise serves to illustrate an interesting point: 
without a finite ^13, equilibrium between and the asymmetries of and cannot be 
quite exact. 



IV. CONCLUSION 

We have given an analytical treatment to the neutrino asymmetry equilibration scenarios 
considered in the numerical studies of DHPPRS 0. In a two-flavour study of the solar 
neutrino oscillation parameters, we demonstrated that the equilibration mechanism is based 
upon a collective adiabatic MSW-like transformation between Ue and (and between z/g 
and I'x), where may be u^, v^, or a linear combination thereof, for synchronised z/g 
and z/g oscillations. The transition temperature is determined, as usual, by Am^^jj and 
by a characteristic momentum p ~ vrT for small initial asymmetries (c.f. (p) ~ 3.15 T for a 
Fermi-Dirac distribution with zero chemical potential). Thus from the sizes of the various 
possible Amg^^'s alone, transition prior to BBN can only be achieved in the case of the LMA 
solution. 

It turns out that the characteristic momentum p does in fact grow with the initial chemical 
potentials of the ensemble, and extremely large |^^|'s can in principle delay the equilibrating 
transition. However, the dependence is weak. For |^^| not exceeding ~ 4, the shift in the 
transition temperature is virtually unobservable. 

A second concern for the LMA solution is the extent of the Ly^/L^^ equilibrium . We 
showed in the present work that, when the mixing parameter is chosen to be maximal 
(as was done in DHPPRS), complete equilibration of these asymmetries can always be 
achieved irrespective of coUisional effects, and for most initial conditions. (The latter can 
in principle alter the adiabaticity of the transition.) For large but not maximal mixing, 
however, the extent of the equilibrium will depend on how one handles collisions on active- 
active neutrino oscillations, and rigorous work on this front is still wanting. Nonetheless, 
heuristic considerations of coUisional damping suggest that full equilibrium before BBN is 
most likely. 

For other solutions of the solar neutrino problem, we demonstrated in a semi- quantitative 
three-flavour analysis that a partial equilibration of L^^ and Lj,^, where Vy is some other linear 
combination of z/^ and z/^, due to a sizeable is possible in principle. Again, the extent of 
this equilibrium hinges on the treatment of scattering processes with the background, and 
is therefore subject to the usual caveats regarding active-active collisions. Investigations on 
this topic are currently underway. 

(Note: See also Ref. for a similar explanation of the DHPPRS results.) 
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